Electronic and Structural Properties of a 4d-Perovskite: 
Cubic Phase of SrZrQ 3 
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First-principles density functional calculations are performed within the local density approxima- 
tion to study the electronic properties of SrZrC>3, an insulating 4d-perovskite, in its high-temperature 
cubic phase, above 1400 K, as well as the generic 3d-perovskite SrTi03, which is also a d°-insulator 
and cubic above 105 K, for comparison reasons. The energy bands, density of states and charge den- 
sity distributions are obtained and a detailed comparison between their band structures is presented. 
The results are discussed also in terms of the existing data in the literature for both oxides. 

PACS numbers: 71.15.Mb, 71.20.-b 



I. INTRODUCTION 

The class of transition metal oxides constitutes a big 
family of interesting materials with extra physical proper- 
ties due to the additional d-electrons they posses. They 
come in variety of crystal structures and exhibit indi- 
vidually several of these phases. They include insula- 
tors, metals, semiconductors and also superconductors. 
Some have delocalized d-bands providing catalytically ac- 
tive surfaces, some have narrow d-bands with emphasized 
electron correlations giving rise to diverse properties like 
high temperature superconductivity, and colossal magne- 
toresistance. They are well known with their ferroelec- 
tric, antiferroelectric and piezoelectric properties. Their 
use in technological application is also diverse, includ- 
ing optical wave guides, laser-host crystals, high tem- 
perature oxygen sensors, surface acoustic wave devices, 
non-volatile memories, dynamic random access memo- 
ries, frequency doublers, piezoelectric actuator materials, 
and high-K capacitors in various applications. 

Strontium titanate, SrTiC>3, is a generic representative 
of transition metal oxides which have perovskite crys- 
talline structure. It has been extensively studied both 
theoretically and experimentally because of its several 
interesting physical and technological properties. It is 
highly insulating at room temperature, and in the form of 
n-type thin films it shows superconductivity at low tem- 
peratures. It is a cubic perovskite at room temperature 
with a tetragonal phase transition at 105 K. Its surfaces 
are very flat and stable both mechanically and chemically 
which makes it best electrode in photocatalysis of wateni, 
and makes it best buffer layer for the growth of gallium 
arsenide on silicon^, and makes it best substrate for the 
growth of high T c cuprate superconductors- . Due to its 
high dielectric constant it is also one of the leading can- 
didates to replace the silica as a gate material in silicon 
technology. 

Another insulating perovskite is the strontium zir- 
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conate, SrZrC>3, with 4d-electrons, which is of interest 
because of its high temperature electronic properties. 
Large single crystals of SrZrC>3 with high perfection can 
be grown with recent techniques and this enables their 
usage as laser-host materials and as substrate materials. 
It was also suggested by Shende et &1& that these mate- 
rials can be used in high-voltage capacitor applications 
because of their high breakdown strengths as well as high 
dielectric constant. In addition to this, both SrTiC>3 and 
SrZr03 are suitable for use in high-temperature applica- 
tions such as fuel cells, steam electrolysis and hydrogen 
gas sensors£*§*£. This is because when these type of tran- 
sition metal oxides are doped with acceptor ions they 
exhibit protonic conduction at high temperatures^. 

Unlike SrTiC>3, at room temperature SrZr03 has an 
orthorhombic phase as revealed by the first structural 
studies which date back to 1960's2ii£. Later the existence 
of two additional phases at high-temperature was pro- 
posed by Carlsson to be both tetragonalii, however, more 
recent studiesi2ii2ii4 on high temperatures have shown 
that SrZrC>3 undergoes three structural phase transitions 
summarized as follows: First, orthorhombic (Pnma) 
to orthorhombic [Cmcm) at 970 K, then to tetragonal 
(I A/ mem) at 1100 K, and then to cubic (Pm3m) at 1400 
K. This compound has a rather high melting temper- 
ature of about 2920 Kl£, consequently it is cubic in a 
wide range of temperature where most of its useful ap- 
plications take place. 

In this work we have made first-principles pseudopo- 
tential calculations of the electronic band structure, den- 
sity of states and charge densities for SrZrC>3 in the cubic 
perovskite phase. In addition to that we have made a ref- 
erence calculation for SrTi03 in order to give a discussion 
by comparison. We have also made some comparisons 
with the related experimental data where available. 



II. CALCULATION METHOD 

We used pseudopotential method based on density 
functional theory in the local density approximation 
(LDA). The self consistent norm conserving pseudopo- 
tentials are generated by using the Troullier-Martins 
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schemed which is included in the fhi98PP packaged. 
Plane waves are used as a basis set for the electronic 
wave functions. In order to solve the Kohn-Sham equa- 
tions^, conjugate gradients minimization method^ is 
employed as implemented by the ABINIT code^S. The 
exchange-correlation effects are taken into account within 
the Perdew-Wang scheme^ as parameterized by Ceperly 
and Alder—. 

Pseudopotentials are generated using the following 
electronic configurations: For Sr 5s electrons are consid- 
ered as the true valence. Moreover the 4s and Ap semicore 
states are added to the valence states. For O only the true 
valence states (2s and 2p) are taken into account, because 
these states are enough to have the correct transferability 
property. For Ti 4s and 3d true valence states plus the 
3s and 3p semicore states are treated as valence states. 
Similarly, for the same group element Zr, 5s and Ad true 
states and additionally the 4s and Ap semicore states are 
considered as valence states. Inclusion of semicore states 
in the case of Sr, Ti and Zr is required in order to get 
the correct electronic properties of these elements in var- 
ious physical systems. In other words, the inclusion of 
these semicore states make the corresponding pseudopo- 
tentials closer to the all-electron potentials. The above 
configuration is found to be the optimized choice for these 
materials. 

All of the calculations involve 5-atom cubic unit cell 
arranged in a perovskite structure. We get a good con- 
vergence for the bulk total energy calculation with the 
choice of cut-off energies at 30 Ha for SrTiC>3 and at 33 
Ha for SrZrC>3 using 4x4x4 Monkhorst-PackSi mesh 
grid. We have found that in the band structure calcu- 
lations 76 k-points are enough to obtain good results 
for both of these transition metal oxides. In the density 
of states calculations, however, the irreducible Brillouin 
zone was sampled with 560 and 455 k-points for SrTiC>3 
and SrZrC>3, respectively. 



III. RESULTS AND DISCUSSION 

The results for structural parameter calculations are 
summarized in Tabic [I] The calculated lattice param- 
eters for SrTiC>3 and SrZrC>3 are both within 0.5% of 
the experimental results. Likewise the calculated bulk 
moduli for SrTiC>3 and SrZrC>3 are found to be 0.3% 
and 0.8% smaller than the effective experimental values, 
respectively. The agreement with the experiments can 
be considered to be very good. In the case of SrTiC>3, 
comparisons of our results with the theoretical work of 
Kimura et ai— and with the calculated values of van 
Bcnthem et ai— suggest that our pseudopotentials are 
as reliable and perform slightly better. To our knowledge 
no first-principles calculation is available for SrZrC>3 to 
compare with. 

For the bulk modulus of SrZrC>3, we have listed in Ta- 
ble H] an extrapolated value of 150 GPa by Ligny and 
Richet— , however, such extrapolation is known to give 



about 15% underestimation for the Zr-compound in a 
series of CaMOg^, and by the same token we expect the 
experimental value to be higher than 150 GPa. 





Lattice Parameter (A) 


Bulk Modulus (GPa) 




Calc. Exp. 


Calc. Exp. 


SrTiOs 


3.878 3.905 


191 183^i 


SrZrOs 


4.095 4.109 — 


171 150 — 



TABLE I: Calculated and experimental values for lattice pa- 
rameter and bulk modulus of SrTi03 and SrZrOs 

Energy band structures and densities of states are 
given in Fig. ^ where the zero of energy is chosen to 
coincide with the top of the valence band. The general 
features of the energy bands are similar for both oxides. 
An overall look at the two band structures shows that 
the lower valence bands are composed of O 2s and Sr Ap 
semicore states grouped together at about -15 eV. Al- 
though the individual bands have similar bandwidths for 
both materials, those of SrTiOs do not overlap and sep- 
arated by a ~1 eV gap, whereas those of SrZrOs having 
stronger interaction with each other, especially around X 
point, causing an overlap of the corresponding densities 
of states, and consequently, the combined bandwidth is 
smaller. 

The upper valence bands have the same trend, i.e., 
wider for SrTi0 3 (4.93 eV) and narrower for SrZr0 3 (4.32 
eV). Top of the valence bands reflect the p electronic 
character mostly due to oxygen-oxygen interaction (down 
to about -4 eV) for both of the transition metal oxides. 
This agrees well with the DV-Xa molecular orbital study 
of Yoshino et al^ '-. Even though whole valence band 
structure is dominated by O 2p states, there are, however, 
a mixture of cr-bands stemming from the pda interactions 
extending throughout the whole valence band width, and 
7r-bands due to pdix interactions which are narrower than 
that. 

Both of SrTiOs and SrZrOs have their valence band 
maxima at the R point. The energy values of the upper- 
most band at T and M points lie slightly lower than its 
value at the point R. In SrTiOs bottom of the valence 
band occurs at R, whereas in SrZrOs it occurs at M point. 
In both cases, however, these energy eigenvalues are close 
to each other at M and R points with a difference, larger 
one being for the SrZrOs, of only 0.25 eV. 

The lower conduction bands are mainly Ti 3d or Zr 
Ad states hybridized with some O 2p electrons giving rise 
to antibonding ir*- and er*-bands. In the case of SrTi03 
the conduction band ordering is as follows: Ti 3c? t2 S (7T*- 
triplet) bands stand alone next to the gap with no overlap 
to the Ti 3d e g states (cr*-doublet) which lie just above, 
and upper parts of which are mixed with the lower exten- 
sion of Sr 5s and Ad t2 g bands. This ordering is slightly 
different for SrZrOs, especially at high energies. Zr Ad t^g 
states are next to the gap and their upper parts are mixed 
with Sr 5s and Ad X,i g states. Zr Ad e g states lie further up 




in the energy region between 8.1 eV and 12.5 eV. Conduc- 
tion band minimum occurs at T point in both materials. 
The lowest conduction band along TX is more dispersed 
in strontium zirconate, and consequently the conduction 
band edge of its density of states is of three dimensional 
nature. Calculated 7r*-conduction band widths are 2.83 
eV for SrTiC-3 and 3.39 eV for SrZr0 3 . 

For SrTi03 the calculated indirect energy band gap be- 
tween r and R points is found to be 1.92 eV and the direct 
band gap is 2.30 eV. Corresponding experimental values 
are 3.25 eV and 3.75 eV^, respectively. The results ob- 
tained for SrZrOs are on the other hand as follows: In- 
direct band gap is 3.37 eV between T and R points. The 
direct band gap being 3.62 eV, is again smaller than the 
experimental value of 5.9 eV— . This disagreement be- 
tween the calculations and experimental values resulting 
in narrower theoretical band gaps is a well known arti- 
fact^ of LDA and does not have any significant effect on 
the rest of the band structure. 

Looking at the density of states pictures one observes 
several structures common to both materials. Most of 
these correspond to singularities in the bands. There are 
three flat bands in the conduction band of both materials. 
Lowest one is Ti (Zr) t2 ff band along XM at about 4 
eV (6 eV) which causes the well defined it* peak with a 
logarithmic van Hove singularity in the density of states 
characteristic of two dimensionality of these bands. The 
next one up is the Ti (Zr) e g band along TX at 4.68 eV 
(8.24 eV) causing a jump discontinuity in the density of 
states. Similar a* shoulder in density of states is caused 
by the third flat band which is again Ti (Zr) e g band 



along MR at about 9 eV (12.5 eV, not shown). 



Valence bands are slightly different in terms of van 
Hove singularities. SrTiC>3 has a very significant flat 
band at -2.87 eV along TXMX which is a non-bonding 
(To band due to O 2p and responsible for the highest peak 
at the center of the valence band. Ti 3d t2 g band at the 
same energy along FX is rather flat, and at about -4.5 
eV along XM and also RX it is quite flat close to the X- 
side. In addition, the band at -4.25 eV along MR looks 
almost flat but the corresponding density of states shows 
three dimensional behavior just slightly. For the SrZr03 
valence band, at about -4 eV, the Zr t2 9 state along XM, 
MR and RX are the only flat bands causing the highest 
peak located at the bottom of the valence band. The 
band corresponding to the oxygen non-bonding state lo- 
cated at the center of SrTi03 valence band density of 
states is dispersed in all directions for SrZr03. 



Charge density plots were obtained for the lowest three 
conduction bands degenerate at the T point for a (001) 
plane containing the transition metal and the four neigh- 
boring oxygens. Fig. [21 shows the composition of these 
bands to be clearly the hybridization between transition 
metal d-orbitals with t2 g symmetry and O 2p orbitals. 
The corresponding tt* bands shown in Fig. ^ arc more 
singled out for SrTiC>3 and consequently the charge is 
more localized as compared to SrZrC>3 whose it* bands 
are not separated from the rest of the conduction bands. 



FIG. 2: Charge density contour plots for tt* bands for (a) SrTiOs and (b) SrZr03 



IV. CONCLUSION 

The structural and electronic properties of two d°- 
insulator metal oxides, SrZrC>3 and SrTiC>3, with cubic 
perovskite structure are studied using an ab-initio pseu- 
dopotential method. Structural parameters are found to 
compare well with the available data in the literature. 
A detailed description of their energy bands are given. 
Corresponding density of states are presented and the 
major structures in them are identified. Charge density 
functions are displayed for the lower conduction bands 
for both oxides. Our results for the electronic properties 



of SrTiOs are shown to agree with other calculations and 
experimental findings whereas those of SrZr03 are com- 
pared only with the estimations of Lee et alm& from their 
optical conductivity spectra, since to our knowledge, no 
theoretical calculations exist in the literature. 
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